Role of symmetry and dimension for pseudogap phenomena 
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The attractive Hubbard model in d = 2 is studied through 
Monte Carlo simulations at intermediate coupling. There is 
a crossover temperature Tx where a pseudogap appears with 
concomitant precursors of Bogoliubov quasiparticles that are 
not local pairs. The pseudogap in A{\iF,uj) occurs in the 
renormalized classical regime when the correlation length is 
larger than the direct ion- dependent thermal de Broglie wave 
length, £^th = hvF (k) /ksT. The ratio TxjTc for the pseudo- 
gap may be made arbitrarily large when the system is close 
to a point where the order parameter has SO (n) symmetry 
with n > 2. This is relevant in the context of SO (5) theories 
of high Tc but has more general applicability. 

71.10.Fd, 71.10.Pm, 71.10.Hf, 74.20.Mn 



Normal state pseudogaps observed in angle-resolved 
photoemission experiments (ARPES) and tunneling 

have become the focus of much of the research on high- 
temperature superconductors. Theoretically, strong- 
coupling models of this phenomenon include repulsive 
doped insulator models H, and attractive models with 
preformed local pairs |^^. At intermediate coupling, 
resonant pair scattering has been invoked |7j. In weak 
coupling, antiferromagnetic or superconducting fluctua- 
tions also lead to pseudogap formation before long range 
order is established [p|-pd|. Finally, arguments about 



the relevance of phase fluctuations |12| do not specify 
whether the small superfluid density comes from strong 
coupling or from thermal fluctuation effects. 

Although the effect of superconducting fluctuations on 
the density of states was studied long ago ||T^, the rel- 
evance of these critical fluctuations on pseudogap phe- 
nomena in high temperature superconductors has been 
questioned mainly because, even in low dimension, the 
critical region should be quite small compared with the 
size of the pseudogap region observed in high temper- 
ature superconductors. Indeed, even for the Kosterlitz- 
Thouless (KT) transition in two dimensions, the range of 
temperatures where superconducting phase fluctuations 
occur is of order ksT'^jEp unless disorder depresses 
[Q . Another question for the original critical-fluctuation 
calculations is that, contrary to density of states pseudo- 
gaps, a pseudogap can appear in the single-particle spec- 
tral weight yl(ki?,Li;) only when the fluctuation- induced 
scattering rate at the Fermi surface ImS {\i.p.,LL} = 0) in- 
creases with decreasing temperature, a behavior that is 



opposite to that predicted by the traditional phase-space 
arguments of Fermi liquid theory. 

The last objection has already been answered |^ by 
showing that in d — 2 the Fermi liquid arguments fail 
when one enters the so-called renormalized classical (RC) 
regime of fluctuations (d = 3 is the upper critical dimen- 
sion H] ) . This RC regime is the one where critical slowing 
down leads to increasing dominance of classical fluctua- 
tions (TiWc < ksT) as temperature is lowered. However, 
the question of the large temperature range where pseu- 
dogaps appear has to be addressed theoretically both at 
the qualitative and at the quantitative levels. On the 
quantitative side, the dependence of the critical region 
on interaction strength and quasi two-dimensionality is 
not obvious, but it has been argued jl^ that specific 
models can give large critical regions even at intermedi- 
ate coupling. On the qualitative side, a factor that may 
considerably enlarge the size of the pseudogap regime is 
the proximity to a point where the order-parameter sym- 
metry is 5*0 (n) with n > 2. Indeed, in d = 2 the critical 
region can become much larger when n > 2 since then 
the Mermin- Wagner theorem implies that the transition 
temperature is pushed down to T = 0. We argue that 
both conditions, namely d = 2 and higher symmetry, are 
generic for high-temperature superconducting materials. 
Indeed, in the underdoped region, where the pseudogap 
is largest, these materials are highly anisotropic (quasi 
two-dimensional) and it has been proposed that the or- 
der parameter may have both antiferromagnetic and su- 
perconducting character corresponding to approximate 
50(5) symmetry @. 

The attractive Hubbard model may be used to illus- 
trate the properties of pseudogaps that appear in such 
situations of approximate high-symmetry in d = 2. For 
our purposes, the important characteristic of this model 
is that the long-wavelength critical behavior is as fol- 
lows for all values of interaction U. At half-filling, there 
is a zero-tcmperature phase transition that breaks the 
finite-temperature SO (3) symmetry JT?! while away from 
half-filling, there is a KT transition at finite temper- 
ature. The corresponding ground state breaks 5*0 (2) 
symmetry. While the details of this model are clearly in- 
appropriate for high-temperature superconductors, it is 
useful to illustrate a number of general points that should 
be applicable to models with transition temperatures 
that are pushed down from their mean-field value by a 
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combination of low dimension and high order-parameter 
symmetry SO{n>2) 

In this paper, we present Monte Carlo simulations for 
the d = 2 Hubbard model at U = —At, a value that is 
slightly on the BCS side of the BCS to Bose-Einstein 
crossover curve {U <U {T™'^^)) Q. We use units where 
nearest- neighbor hopping iat = 1, lattice spacing is unity, 
fi = 1 and fcs = 1. Previous numerical work charted 
the phase diagram |Q . They have also investigated the 
pseudogap phenomenon mostly in strong coupling where, 
we stress, the Physics is different from the case discussed 
here [|I|,|o|,|l) . On the weak-coupling side of the BCS to 
Bose-Einstein crossover, there have been numerical stud- 
ies of KT superconductivity as well as several discus- 
sions of pseudogap phenomena in the spin properties and 
in the total density of states at the Fermi level |^,^,^ . 
The only study of A (k, lo) was restricted to regions far 
from the SO (3) symmetric point ||2l[| . 

Here we establish a dynamical connection between the 
appearance of the RC regime in the pairing collective 
modes and pseudogap formation in single-particle quan- 
tities. In particular, we show that a) Close to a high 
symmetry point, pseudogaps can appear at a crossover 
temperature Tx that scales with the mean-field transi- 
tion temperature while the real transition may occur at 
much lower temperature, Tc, leading to a wide tempera- 
ture range for the pseudogap. b) At the crossover temper- 
ature, one enters the RC regime where the characteristic 
frequency for fluctuations becomes smaller than the tem- 
perature, c) Pseudogap in weak-to-intermediate coupling 
do not require resonance in the two-particle correlations 
p^ . d) To have a pseudogap in A (k^,!^) for a given wave 
vector, it is not enough to have the collective mode (two- 
particle) correlation length satisfy ^ > 1. It is necessary 
that ^ becomes larger than the single-particle thermal 
de Broglie wavelength £^th = vp (k) /T. This implies in 
particular that even for an isotropic interaction, as tem- 
perature decreases a pseudogap opens first near the zone 
edge, where vp is small, and it opens last along the zone 
diagonal where vp is largest. This anisotropy would be 
amplified for an anisotropic interaction of d— wave type 
Ip^ . For U — —4, the condition ^ > S^th is realized near 
the zone edge for ^ not so large. Analytical arguments 
for the above results have appeared elsewhere ||Jl^,^ . 

Let us first recall a few results at half-filling, (n) = 1, 
where the chemical potential fi vanishes. There the 



canonical transformation Cii — > ( — 1) 



maps the 



attractive model onto the repulsive one at the same fill- 
ing. The q = 0, s— wave superconducting fluctuations 
and the Q — (tt, tt) charge fluctuations arc mapped onto 
the three components of antiferromagnetic spin fluctua- 
tions of the repulsive model and hence they are degen- 
erate. Because of this degeneracy, the order parameter 
at half-filling has 5*0(3) symmetry hence, by the 
Mermin- Wagner theorem, in two dimensions the phase 
transition is at = 0. Results for the attractive model 
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FIG. 1. Pair (filled symbols) and charge (open symbols) 
structure factors for U = —4 and various system sizes, (a) 
Filling (n) = 0.95. (b) Filling (n) = 0.8. 



are easily extracted from simulations of the canonically 
equivalent repulsive model. The pair structure factor 
5a = (AtA -t- AAt) with A = E,=i cnca and the 
Q=(7r,7r) charge structure factor 5c = (pqP-q) are 
identical, showing an increase as temperature decreases 
and then size-dependent saturation. The sudden rise of 
S'a as temperature decreases indicates a crossover to a 
RC regime with a concomitant opening of the pseudogap 
|p|,p|j2^ in A {}ip,uj). The crossover temperature is a size- 
able fraction of the mean-field transition temperature. 

Slightly away from half-filling, the SO (3) symmetry 
is formally broken by the chemical potential since it 
couples only to the charge part of the triplet. However, 
in the regime where the temperature is larger than the 
symmetry-breaking field, the symmetry is approximately 
satisfied. For filhng (n) = 0.95, we have T > 
for the whole range of temperature shown in Fig. 1(a) 
[jjL ^ —0.07 at T = 0.125). On this figure, superconduct- 
ing (S'a, filled symbols) and charge (S'c, open symbols) 
structure factors are of comparable size when one en- 
ters the RC regime, showing that we have approximate 
SO (3) symmetry for various sizes L x L. The beginning 
of the RC regime, that occurs at a temperature nearly 
identical to the crossover temperature j|] Tx identified at 
half-filling, is signaled by the increase in the magnitude 
of correlations. Eventually, the concomitant increase of 
^ leads to the size-dependence apparent at lower tem- 
perature in Fig. 1(a). The plot in Fig. 1(a) resembles the 
result at half-filling. [|^,| The near equality of supercon- 
ducting and charge fiuctuations at the crossover to the 
RC regime should be contrasted with the case (n) — 0.8 
in Fig. 1(b) where the charge fluctuations show basically 
no critical behavior when the superconducting correla- 
tions begin to do so. Hence, at this flUing, there is little 
SO (3) symmetry left at Tx ■ This is expected since the 
symmetry-breaking field |/x(Tx)| — |/x(0.25)| = 0.26 is 
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FIG. 2. Even part of the pair spectral weight Xa('^)/'^ ob- 
tained from analytic continuation of imaginary time Monte 
Carlo data for U = -4, Ar = 1/10, (n) = 0.95 and size 8x8, 
except for symbols that are for 6 x 6 at T = 1/5. About 10^ 
sweeps were done in each case. 



comparable to Tx- One basically enters directly into the 
RC regime of a SO (2) KT transition In this regime 
dTc {n) /dn > 0. 

Let us go back to the filling (n) = 0.95. At this fiUing, it 
has been estimated that Tc < 0.1 and dTc (n) /dn < 
0. We already showed that at this filling, one has ap- 
proximate SO (3) symmetry. One can confirm that the 
increase of and Sc at Tx comes from RC fiuctua- 
tions by calculating the spectral weight for superconduct- 
ing fluctuations x'a (imaginary part of the T — matrix) , 
Xa i^) = Jdt^ l (t) ' ) e*"*- The even part H of 
Xa ('^) /'^ obtained from a Monte Carlo calculation 
of the imaginary-time quantity (A (t) A''^ -|- A^ (r) A) 
followed by Maximum Entropy inversion p9| ]. The even 
part of x'a i'^) /"^ is plotted on Fig. 2 for an 8 x 8 system 
and various temperatures. The symbols are for a 6 x 6 sys- 
tem at T — 1/5. The maximum is at zero frequency, as 
for an overdamped mode Xa^ (Q;^) °^ S, '^ +9^ ^ iuj/tuo. 
The characteristic frequency, given by the half-width at 
half maximum, is ujc = wq/^^ with ujq a microscopic re- 
laxation rate. There is a marked narrowing of the width 
in frequency as temperature decreases. One enters the 
RC regime when — 1/4 ~ Tx, a temperature larger 
than Tc(< 0.1). At T = 1/5, the correlation length is 
becoming comparable with system size since the 6x6 
system gives a result that differs from 8x8. 

The effect of RC collective fluctuations on single- 
particle quantities is illustrated in Fig. 3 (a) to (c) that 
show density plots of the single-particle spectral weight 
A (k,u;) for an 8 x 8 system at, respectively, T — 1/3, 1/4, 
and 1/5. When temperature reaches T = 1/4 one no- 
tices that a minimum (pseudogap) around k = (0,7r), 
Lu ^ develops along with two maxima away from 
Lu = 0. The latter maxima are precursors of the Bo- 
goliubov quasiparticles of the ordered state ||^ |Q . The 
pseudogap becomes deeper and deeper as temperature 
decreases, the distance between maxima remaining about 
constant, as observed in high temperature superconduc- 
tors 1^. The condition for the appearance of a pseu- 
dogap in A (kjw) is not only that we should be in the 
RC regime and in low dimension but also that ^ /^th 
should be large This illustrated by the fact that 
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FIG. 3. Density plot of the single-particle spectral weight 
A{]i,uj) from about 10^ Monte Carlo sweeps for U = —4, 
At = 1/10, (n) = 0.95 and size 8x8. Top to bottom, 
T = 1/3,1/4,1/5. The dilation by ^(2) of the axis from 
(tt, tt) to (0, 0) allows a comparison of Fermi velocities. 



the pseudogap at k = (7r/2, 7r/2) , where the Fermi veloc- 
ity is larger, is not opened yet at T = 1/4 despite the 
fact that at T = 1/5 the pseudogap is of comparable 
size both around k = (0,7r) and at k = (7r/2, 7r/2), in 
concordance with the equality of the gaps at these two 
points in the zero-temperature spin-density wave state. 
From the slopes in Figs. 3 (a) to (c), vp is clearly larger at 
(7r/2, 7r/2) meaning that the condition ^ > £,th = vp/T 
is harder to satisfy at this wave vector. Numerical es- 
timates show that ^ is nearly isotropic by contrast with 
vp- These estimates are consistent with the appearance 
of the pseudogap in A (k,w) when ^ ~ ^th- 

As in any numerical simulation, finite-size effects 
should be considered carefully. There are two impor- 
tant intrinsic lengths in this problem, namely S,th and ^. 
When L ^ £^th, the system acts as if it was basically 
in the quantum zero-temperature limit of a finite sys- 
tem. We have checked that at T = 1/8, A(k, w) shows 
real gaps, instead of pseudogaps, that appear at progres- 
sively higher temperature in systems of smaller size. For 
T — 1/3, 1/4 on the other hand, estimates oi £,th and of ^ 
as well as calculations for 6x6 systems and a few 10 x 10 
systems suggest that our numerical results for A (k,uj) on 
8x8 systems are free of appreciable size effects, i.e. L 
> S,, ith- This can also be checked by the size dependence 
of the results in Fig. 1(a) and Fig. 2. While size effects be- 
come important in x'a ('^) 1^ when ^ exceeds L, as long 
as £,th < L then it is possible to see thermally induced 
pseudogap effects in A (k, to) even if ^ > L [||. 

In summary, the qualitative phase diagram for the at- 
tractive Hubbard model in d = 2 sketched in Fig. 4 shows 
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FIG. 4. Schematic crossover diagram for the d = 2 attrac- 
tive Hubbard model in the weak-coupling regime. The shaded 
area above the KT phase, up to the dotted, line is the RC 
pseudogap regime. The KT critical regime is hatched. 



that near a point with high order parameter symmetry, 
the transition temperature decreases while the pseudogap 
temperature increases along with the mean-field transi- 
tion temperature and the zero-temperature gap. In the 
region where dTc/dn < 0, the crossover to KT critical 
behavior occurs in the RC pseudogap regime when T 
is less than the symmetry breaking field. In our sim- 
ulations we would need larger system size to reach the 
KT critical regime. Contrary to the scenario of Ref. 
0, in our case dimension and symmetry contribute to 
create a wide pseudogap region, there is no critical cou- 
pling strength, and furthermore one enters the RC regime 
without sharp resonance in x'a i^) 1^- Also, the precur- 
sors of Bogoliubov quasiparticles in A (k, iS) occur under 
conditions very different from those for strong-coupling 
Cooper pairs that are local and do not need low dimen- 
sion or large ^/^t/i. Comparisons with non-perturbative 
many-body calculations should appear elsewhere js^ . 

In high Tf. superconductors the competition is between 
antiferromagnetism and superconductivity. Recent time- 
domain transmission spectroscopy experiments [ |3l[ | sug- 
gest that the RC regime for the KT transition (hatched 
region in Fig. 4) has been observed. Close enough to the 
transition there is dimensional crossover to c? = 3. For an- 
tiferromagnetic fluctuations, there are suggestions from 
NMR of a RC regime , but there is no definite proof. 
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